Where are the missing baryons in clusters? 
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Observations of clusters of galaxies suggest that they contain significantly fewer baryons (gas 
plus stars) than the cosmic baryon fraction. This 'missing baryon' puzzle is especially surprising 
for the most massive clusters which are expected to be representative of the cosmic matter content 
of the universe (baryons and dark matter). Here we show that the baryons may not actually be 
missing from clusters, but rather are extended to larger radii than typically observed. The baryon 
deficiency is typically observed in the central regions of clusters (~0.5 the virial radius). However, 
the observed gas-density profile is significantly shallower than the mass-density profile, implying 
that the gas is more extended than the mass and that the gas fraction increases with radius. We 
use the observed density profiles of gas and mass in clusters to extrapolate the measured baryon 
fraction as a function of radius and as a function of cluster mass. We find that the baryon fraction 
reaches the cosmic value near the virial radius for all groups and clusters above ~ 5 x lO^'^/ifj^MQ. 
This suggests that the baryons are not missing, they are simply located in cluster outskirts. Heating 
processes (shock-heating of the intracluster gas, plus supernovae and AGN feedback) that cause the 
gas to expand are likely explanations for these results. Upcoming observations should be able to 
detect these baryons. 
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I. INTRODUCTION 

Clusters of galaxies, the largest virialized systems in 
the universe, are powerful tools in constraining cosmol- 
ogy and tracing the large-scale structure of the universe 
[1|-1J, and references therein] . The large mass of clusters 
(~ 10^'* to 10^^h^2-^Q) implies that their contents — dark 
and baryonic matter — have been accreted from very large 
regions of '--^lO comoving Mpc, and therefore should be 
representative of the mean matter content of the uni- 
verse. The strong gravitational potential of clusters also 
implies that baryons cannot easily escape from these sys- 
tems. Therefore, clusters are expected to retain the cos- 
mic baryon fraction, i.e., the relative fraction of baryons 
to total matter. This basic assumption of a representa- 
tive baryon fraction in clusters was used in 1993 Q to 
suggest that the mass-density of the universe must be 
low, since the observed baryon fraction in clusters was 
considerably larger than expected for a critical density 
universe. Most of the baryons in clusters reside in the X- 
ray emitting hot intracluster gas, which approximately 
traces the cluster gravitational potential dominated by 
dark matter. The rest of the baryons are in the luminous 
galaxies and in isolated stars that comprise the small 
amount of faint diffuse intracluster light (ICL). 

A puzzle has been raised, however, over the last few 
years: Detailed X-ray observations from Chandra, XMM- 
Newton and others suggest that the cluster baryon frac- 
tion (gas plus stars relative to total mass) is considerably 
lower than the cosmic value |6l-(ll|. The cosmic baryon 
fraction is well determined both from Big-Bang nucle- 
osynthesis [ij, [iJl and from observations of the Cos- 
mic Microwave Background to be /;, = 0.1675 ± 0.006 
[WMAP7: (H]. The cluster gas fraction has been re- 
ported by the above observations to be only 60 — 80% 
of the cosmic value, with stars contributing only a small 
(^10%) additional amount of baryons. This baryon dis- 



crepancy is observed to increase with decreasing cluster 
mass [i3, [ill- This raises the questions: Where are the 
missing baryons? Why are they 'missing' ? 

Attempted explanations for the missing baryons in 
clusters range from preheating or other energy inputs 
which expel gas from the system |15l - [2ll and references 
therein] , to the suggestion of additional baryonic compo- 
nents not yet detected [e.g., cool gas, faint stars la. l2a| . 

In this paper we investigate the possibility that the 
'missing baryons' are not missing at all, but are rather 
located in the outskirts of clusters where few detailed 
observations have yet been made. The 'missing baryons' 
problem is typically observed within the central regions 
of clusters, generally within a radius of i?5oo (where the 
enclosed mass-density is 500 times the critical density). 
This radius is ^0.5 of the virial radius of the cluster 
(where the enclosed density is ^100 times the critical 
density, for the current LCDM cosmology [23, 2^ ) ■ Thus 
for a virial radius of ^1.5 Mpc, the typical missing baryon 
problem is observed only at ^0.75 Mpc from the cluster 
center. 

Observations show that the gas density profile in the 
outer parts of clusters decreases with radius slower than 
the mass profile in these regions. Using gravitational 
lensing, the latter has been observed out to large radii 0, 

t [39{ and is consistent with the expected NFW profile 
. While the cluster mass-density declines with radius 



approximately as 



-2.6 



in these outer regions, the gas- 



density declines only as r~^-^. This implies that the gas 
is more extended than the total mass, and therefore the 
gas fraction increases with radius beyond the observed 
radius of i?5oo- A shallow slope of the gas profile (as 
compared with the mass profile) is indeed expected if gas 
heating occurs in the clusters (e.g., from shock- heating 
of the gas, supernovae, and AGNs; these processes are 
known to exist). The heating expands the gas relative 
to the dark matter potential and spreads it out to larger 



radii, with a shallower slope, as observed. 

We use the observed slopes of the gas-density and 
mass-density profiles in the outer regions of clusters to 
extrapolate the observed gas fraction from R^qq to larger 
radii, up to the virial radius {Rmr — Rioo [23,[2J|). We 
add the observed stellar fraction to the extrapolated gas 
fraction to find the baryon fraction at large radii. We 
perform this extrapolation as a function of cluster mass, 
from groups to rich clusters. Our analysis is based en- 
tirely on observations. 

We find that there is no 'missing baryon' problem in 
rich clusters when the data is extrapolated to near the 
virial radius; at that radius, the baryon fraction is consis- 
tent with the cosmic value. Most of the 'missing baryons' 
are therefore expected to be in the outskirts of clusters, 
between i?5oo and Rvir- This result can be tested with 
upcoming observations of the Sunyaev-Zeldovich effect in 
clusters (e.g., SPT ^■, ACT [26|]) as well as with more 
sensitive X-ray observations. 

Observations have shown that the 'missing baryon' 
problem at i?5oo becomes significantly more severe for 
lower mass clusters and groups of galaxies than for rich 
clusters; the observed gas fraction decreases considerably 
with decreasing cluster mass. This too would be expected 
if the heating processes expand the gas: the lower gravi- 
tational potential of the smaller systems will not be able 
to hold on to their gas as well as the higher mass clusters. 
The gas density profile in small groups is indeed observed 
to be shallower than the gas density profile in massive 
clusters, suggesting that the gas in low mass systems is 
more spread out. We extrapolate the observed baryon 
fraction as described above for clusters as a function of 
their mass — from groups to rich clusters. We find that 
for this entire mass range the baryon fraction within R^r 
is flat and is consistent with the cosmic value. 

In the next section we describe the observations and 
analysis. A discussion of the results and our conclusions 
follow. We use a LCDM cosmology with h=0.72 and 
a«=0.258. 



radius (and occasionally beyond). The observed stellar 
mass fraction, including the small contribution from the 
ICL, is added to the gas fraction to yield the total baryon 
fraction. The baryon fraction is then investigated as a 
function of radius, from -R500 to Rmr, and as a function 
of mass, from groups to rich clusters. 



A. Gas Fraction at R500 

Although some observations extend to i?200i the gas 
fraction has been accurately measured for a sufficiently 
large sample of nearby clusters only out to i?5oo • 

We use X-ray observations of the gas fractions for 
39 nearby clusters observed with Chandra and XMM- 
Newton |8l-[ig|. These authors use similar methods of 
data reduction and analysis. We use the compilation 
by Giodini et al. [il| of these groups and clusters above 
the mass of M500 = 3 x 10^^h^2^lQ (where M500 is 
the mass within i?5oo)- After conversion to a common 
cosmology (LCDM with Q^ = 0.258 and Hq = 72 
km/s/Mpc), the three samples have been binned into 
four logarithmically spaced mass bins, from groups to 
rich clusters [ll|. (We do not include the lowest mass 
bin at 10^^h^2 -^0 which contains only two groups with 
large error bars.) Our cluster sample has a mass range 
of 3 X W^^Hj^Mq < M500 < lO^'^/ifa^MQ and a redshift 
range of 0.012 < z < 0.23. The mean observed gas frac- 
tion for each mass bin is listed in Table 1 and shown in 
Figure 1 . The error on the mean is the rms standard de- 
viation divided by y/N — 1. The horizontal bars are the 
mass ranges for the bins. Also presented in Figure 1 is 
the cosmic baryon fraction observed by the WMAP7 mi- 
crowave background measurements, ft = 0.1675 ± 0.006 
[14| . One can see that the cluster gas fraction at -R500 
is significantly lower than the cosmic baryon fraction. 
The gas fraction decreases significantly from rich to poor 
clusters; while rich clusters contain about 12% gas within 
R500, the gas fraction in poor clusters and groups is only 
-6-7%. 



II. OBSERVATIONS, ANALYSIS AND RESULTS 

The gas fraction in clusters of galaxies has been mea- 
sured for a relatively large sample of groups and clus- 
ters within R50Q. The total baryon fraction is obtained 
by adding the stellar mass fraction observed for these 
systems within the same radius. This baryon fraction 
is systematically lower than the cosmic value measured 
by WMAP7; the discrepancy increases with decreasing 
cluster mass [10, [ill and references therein] . Here we in- 
vestigate the possibility that the missing baryons are not 
actually missing from clusters but are rather spread out 
to larger radii, beyond i?5oo- We investigate this scenario 
by extrapolating the observed gas fraction to larger radii, 
from i?5oo up to the virial radius {Rmr — Rioo)^ using the 
mean observed gas and mass density profiles in clusters; 
these density profiles have been measured up to the virial 



B. Stellar Fraction 

The galactic stellar fraction has been measured in 
nearby clusters using multiband optical and infrared sur- 
veys combined with stellar population models. We use 
the results obtained from the COSMOS survey [ii| and 
2MASS survey ^8] for nearby clusters. The combined 
COSMOS and 2MASS sample covers our entire cluster 
mass range, 3 x 10^'^ — lO^^ft-fj^Af©. We bin the observed 
stellar fraction into the same four logarithmic mass bins 
as the gas fraction sample. The mean stellar fraction de- 
clines with cluster mass as JVf-o.37±o.o4 [gee[ll|, Figure 5], 
exhibiting an opposite trend to that of the gas fraction. 

To determine the contribution of the diffuse intraclus- 
ter light (ICL) to the stellar fraction, we use the obser- 
vations by Zibetti \2^ which use the Sloan Digital Sky 



TABLE I: Cluster gas, stellar, and total baryon fractions 
within -R500 



Bin 


< M500 > 


rgas 
J 500 


rstars + ICL 
J 500 


/500 


1 
2 
3 

4 


5.1 X 10" 

1.2 X 10^* 

3.0 X 10^* 

7.1 X 10^* 


0.068 ±0.005 
0.080 ±0.003 
0.103 ±0.008 
0.123 ±0.007 


0.050 ± 0.002 
0.040 ± 0.004 
0.023 ± 0.002 
0.021 ±0.002 


0.118 ±0.005 
0.120 ±0.005 
0.126 ±0.009 
0.143 ±0.007 



Gas, stellar (including ICL), and baryon fraction of clusters 

within R500 for four cluster mass bins: Bin 1 
(0.3 - 0.7 X lO"/if2^M0), Bin 2 (0.7 - 1.7 x lO^^'/ifa^M©), 

Bin 3 (1.7 - 4.2 x lO^/ifj^Me) and Bin 4 

(4.2 - 10 X lO^/ifj^Mo). /f^are averages from Chandra 

and XMM observationslffij. /''*"'■'' and the 10% ICL 

contribution are from [lll.[2a. [29l|. The error bars are l-u 

errors on the mean. 



Survey data to reach to unprecedented cluster-centric dis- 
tances and depth for a large number of stacked clusters 
of various masses. They find the ICL is centrally concen- 
trated, and that on average the ICL contributes ~ 10% 
of the stellar light within the central 500 kpc for all clus- 
ter masses. We add this 10% contribution to the galactic 
stellar fraction discussed above for all clusters. We note 
that the ICL contribution may decline to less than 10% 
when extending to larger cluster radii; but since the ICL 
is a very small contribution to the total baryon fraction, 
this effect has negligible consequences (see Discussion). 
The total stellar fraction for the four mass bins is sum- 
marized in Table 1. It is added to the gas fraction to 
obtain the total average baryon fraction for each mass 
range. The baryon fraction within -R500 is listed in Table 
1 and plotted as a function of mass in Figure 1. The 
deficiency of baryons within R500 relative to the cosmic 
value is clearly seen in Figure 1; the deficiency becomes 
more severe for lower mass clusters. 



C. Gas and Mass Density Profiles 

Extrapolating the observed gas fraction to larger radii 
beyond i?5oo requires the knowledge of the observed gas 
and mass density profiles in these regions. The gas den- 
sity profile has been measured well in the outer parts 
of clusters (i?5oo ~ Rvir) using X-ray observations of 
nearby clusters with ROSAT, Chandra, XMM-Newton, 
and Suzaku. The observed gas profile in the outer re- 
gions fits well to a beta-model with a density slope of 



3/3„ 



agas (where pga 



ex r 



'). We use observa- 



tions with small uncertainties on the gas density slope at 
large radii and which cover our entire cluster mass range 
[a, 31-34]. 

The averaged observed /3- values of the gas density slope 
are presented as a function of cluster temperature in Fig- 



ure 2. The data include measurements for 51 clusters as 
well as average slopes for stacked samples of hundreds 
of optical clusters; the slopes are binned in temperature 
as presented in Figure 2. The error bars on the bin- 
averaged gas slopes are the standard deviation divided 
by y^N- 1. The ROSAT observations of 39 z < 0.25 
clusters by Vikhlinin et al. 31] cover the full outer re- 
gions of clusters, from O.3-R500 to 1.5i?i8o {^ Rvir)- The 
best-fit slope to these outer regions is determined for a 
wide range of cluster temperature, from 2 keV to over 
10 keV. The Chandra observations of 10 massive clus- 
ters by Vikhlinin et al. [g provide /3-fits to gas density 
slopes near i?5oo, nicely consistent with the trend shown 
by the slopes of the previous sample [3i| . For the lowest 
mass bin (M500 « 5 x 10^^h^2^o) "^^ use the observed 
density slope from ROSAT by Dai et al. [32, their rich- 
ness class 1] , who obtain integrated X-ray gas profile for 
stacked samples of hundreds of low-mass optical clusters 
out to Ryir- We also include Bautz et al. 33] who mea- 
sure the gas profile of Abell 1745 to i?2oo using Suzaku 
observations. 

The observed beta slopes presented in Figure 2 are 
consistent with each other, and show a shallower slope 
for lower mass systems than for massive clusters. For our 
four mass bins (Table 1, Figure 1) we find the following 



mean gas density slopes: 



•^gas 



= 3 X Pgas = 1.8 ± 0.2 



for Bin 1, 1.9 ± 0.07 for Bin 2, 2.1 ± 0.02 for Bin 3, and 
2.3 ± 0.02 for Bin 4. The error bars are the standard 
deviation of the average /^-values in each bin divided by 
\/N — 1. When comparing a slope at a given temperature 
bin in Figure 2 to a given mass bin in Figure 1, we use 
the observed M500 — T relation 8]; the results are not 
sensitive to the exact conversion because of the slowly 
varying Pgas^T) relation. 

The final piece required for extrapolating the gas (and 
baryon) fraction to large radii is the total mass density 
profile. This has been observed to follow the NFW pro- 
file [2^ to large radii, by using weak lensing observations 
with the Sloan Digital Sky Survey and other observa- 
tions [3, [3^, [33, and references therein]. We use the 
average observed value of the concentration parameter 
C200 — 5 for our mass range [3^. The NFW profile has a 
mass density slope of am = 2.6 in the radius range from 
i?50o to i?200j and am = 2.7 from i?200 to Rmr (where 
Pm oc r~"™). These slopes are considerably steeper than 
the corresponding gas density slopes, thus yielding an 
increasing gas fraction with radius in cluster outskirts. 



D. Extrapolation and Results 

Using the observed gas density and mass density slopes 
at large radii, we extrapolate the observed gas fraction 
from -R500 as a function of radius up to the virial radius. 
The gas fraction increases with radius as: 






r> R 



500- 



The extrapolated gas fraction from i?5oo to i?200 and to 
Ryir is presented as a function of cluster mass in Figure 
3a. Because the gas density slope is shallower in groups 
than in rich clusters, the increasing trend of gas fraction 
with mass becomes weaker at the outer radii. 

The baryon fraction is presented for the different 
radii — R500, R200, and Ryir — as a function of mass in 
Figure 3b; the baryon fraction within these radii is the 
sum of the gas fraction (Figure 3a) and the stellar mass 
fraction discussed above, including the 10% ICL. We as- 
sume that this fraction remains constant with increas- 
ing radius; the main results do not change significantly 
if this assumption is changed because of the relatively 
small contribution of the stellar component (see Discus- 
sion). The error bar on the extrapolated gas fraction 
is the propagated errors of the gas fraction at -R500 and 
the gas slope used for extrapolation. The error bar on 
the baryon fraction is the combined errors of the gas and 
stellar fractions. 

The results in Figure 3b show that the baryon fraction 
flattens considerably as a function of cluster mass when 
extrapolated to larger radii; this is due to the combined 
effect of the shallower gas density profile in groups, which 
results in more gas in their outskirts, plus the larger ob- 
served stellar mass fraction in groups, which adds more 
baryons in the smaller systems. In fact, at the virial ra- 
dius we find that the baryon fraction is essentially fiat 
from groups to rich clusters, at a level consistent with 
the cosmic baryon fraction. This suggests that there are 
no 'missing baryons': most of the 'missing baryons' are 
likely located in the outskirts of clusters, extending to 
nearly the virial radius. 

The extrapolated baryon fraction is presented as a 
function of radius for two of our mass bins (Bins 2 and 4) 
in Figure 4; the results show the slow but steady increase 
in the baryon fraction with radius. 



III. DISCUSSION 

The gas fraction at i?5oo is often cited to claim that 
clusters contain significantly fewer baryons than the uni- 
versal baryon fraction and therefore exhibit a 'missing 
baryon' problem. Here we show, based purely on ob- 
servational results, that the gas (and baryon) fraction in- 
creases substantially with radius beyond i?5oo; the radius 
-R500 cannot be used as a representative radius for com- 
parisons with the global value. Using the observed gas 
and mass density profiles, we extrapolate the observed 
baryon fraction (gas and stars) as a function of radius 
from -R500 to the virial radius. Since the gas-density is 
observed to decline more slowly with radius than the total 
mass, we find that the average baryon fraction increases 
with radius for clusters of all masses; it reaches the cos- 
mic baryon fraction near the virial radius (i?ioo)- This 
suggests that baryons are not missing in clusters, they 
are simply located in cluster outskirts. 

Recent observations of the Sunyaev-Zeldovich (SZ) ef- 



fect in 15 massive clusters using the South Pole Telescope 
|27| measure the gas density pressure profile in these clus- 
ters as a function of radius up to the virial radius (and 
beyond) ; they are well fit with a beta-model in the outer 
regions. The detection of the gas to these large radii, 
and their observed beta-model slope (when corrected for 
the temperature profile), are nicely consistent with our 
results and the conclusion that the baryons are out there 
between i?5oo and Ryir- Similarly, George et al. [36| use 
Suzaku X-ray observations to trace the gas density pro- 
file in the massive cluster PKS0745-191 up to the virial 
radius, observing the baryons at the cluster outskirts and 
measuring a shallow gas-density profile in these outer re- 
gions. 

The gas density profile is observed to be even shallower 
in low-mass clusters than in massive clusters. This is con- 
sistent with the gas being heated via shocks and feedback 
from supernovae and AGN. This feedback will be more 
significant in low-mass systems when compared to the 
binding energy of the gas. If star formation is more effi- 
cient in groups than in clusters [e.g., [2^, this will further 
increase the gas entropy in these systems, because the 
star formation removes the lowest-entropy gas from the 
intracluster medium, leaving behind gas with higher av- 
erage entropy [37|. A higher stellar fraction also implies 
more supernovae/AGN activity per unit mass. This ex- 
plains why, at -R500, groups have a lower baryon fraction 
than clusters. However, by the virial radius, we are able 
to account for all the baryons expected from the cosmic 
value. 

A few comments on uncertainties are in order. First, 
the gas slopes we use from ROSAT observations iSlI, Is^] 
are the best-fit slopes in the large radial range from i?5oo 
(^ Q.bRyir) to ^ Rvir- This is the average slope in these 
outer cluster regions which are relevant to our extrapo- 
lation. If we instead use the best measured gas-density 
slope at i?5oo [e.g., Q and then steepen it slowly with 
radius similar to the steepening trend observed for the 
NFW profile (which steepens by 10% from i?5oo to Rvir), 
we find consistent results within the error-bars with the 
ones presented above. 

The stellar fraction has been observed for a large num- 
ber of clusters by Gonzalez et al. [30| , who report a some- 
what higher stellar fraction for low-mass clusters and a 
slightly lower stellar fraction for high-mass clusters than 
the observations by Giodini et al. 11 1 and Lin et al. !28| . 
The higher stellar fraction in groups [30] is likely due to 
a selection bias towards systems with dominant Bright- 
est Cluster Galaxies in the small groups. Their observed 
trend of stellar fraction with mass is therefore somewhat 
steeper, log/"*"'^'' = (7.57±0.08)- (0.64 ±0.13) log M500. 
Using this stellar fraction in our analysis does not change 
our result significantly; the baryon fraction at the virial 
radius decreases by « 5% for the most massive bin, and 
increases by « 7% for the lowest mass bin, but all are 
well within our 1-a error bars. 

We use the observed 10% contribution of the ICL to 
the stellar fraction [29]. If the ICL fraction changes at 



the outer radii, it will not affect our results since the ICL 
contribution is only ^0.2-0.4% of total mass. 

Similarly, we assume in our analysis that the stellar 
fraction does not change when extrapolated from -R500 
to Rvir , since no observations of the stellar fraction have 
been made beyond R500 . If the stellar mass fraction de- 
creases somewhat with radius, then the baryon fraction 
will slightly decrease. If we assume, for example, a 20% 
decrease in stellar fraction from i?5oo to Rvir , we find that 
the baryon fraction is lowered by only 2% for the massive 
clusters and 8% for the low- mass groups (at Rmr)] these 
values are within l-cr of our baryon fraction results even 
for the lowest mass groups. 

For the mass density profile, we use the mean observed 
value of the concentration parameter of the NFW profile, 
C200 = 5, as observed from weak lensing (38] for our mass 
range. As Mandelbaum et al. [3g| discuss in their paper, 
this value is slightly lower than results from simulations 
and some previous studies. If C200 is larger than 5, the 
mass profile will be more concentrated, i.e., fall off even 
steeper with respect to the gas profile at the outer radii. 
This would cause the gas fraction to increase even more 
with radius. The effect is small, however, and a change 
of C200 to 7 induces a change in baryon fraction at the 
virial radius of <5%. Earlier in this Section we discussed 
a scenario in which the gas density declines at the same 
trend as the NFW profile; in that scenario, the effects 
will nearly cancel out. 

The observations presented above are consistent with 
an energy input in clusters that heats the intracluster 
gas and expands it to larger radii, where the missing 
baryons may be found. Hydrodynamic simulations of 
cluster formation which include the relevant physics are 
in qualitative good agreement (within ~10%) with the 
picture presented here. Using simulations without cool- 
ing, star formation, or feedback, the baryon fraction is 
roughly constant, at 90% of the cosmic value, from i?5oo 
to i?200 t4^j- When these processes are included (with 
feedback coming from AGN), then the baryon fraction 
is found to be lower at -Rnnn, but instead of being fiat 
it increases with radius [4l|, |42]. Similar results, based 
on cluster energetics, are found in the models of Bode et 
al. |2l| . Further improvements in models and simulations 
(including additional physical processes, sources of non- 
thermal pressure and possible non-equipartition effects) 
will be needed to provide precise comparisons with the 
observations. 



IV. CONCLUSIONS 

We investigate the 'missing baryon' problem in clus- 
ters of galaxies. Observations show that the baryon frac- 



tion (gas plus stars) measured within a radius of i?5oo 
in groups and clusters is significantly below the cosmic 
value; this baryon discrepancy increases with decreasing 
cluster mass. This gives rise to the puzzle: Where are 
the missing baryons? Why are they missing? We show 
that the baryons may not be missing at all, but rather 
are spread out to larger radii, beyond i?5oo, and can be 
found in the cluster outskirts. Here we show that the ra- 
dius i?5oo, which is typically used in observations, cannot 
serve as representative for the baryon fraction in groups 
and clusters. Based entirely on observations, we investi- 
gate the dependence of the baryon fraction on radius for 
clusters of different masses, from groups to rich clusters. 
We use the mean observed gas and mass density profiles 
in clusters to extrapolate the observed gas fraction as a 
function of radius from i?5oo to the virial radius. Since 
the gas-density profile is significantly shallower than the 
mass-density profile, the gas-fraction increases with ra- 
dius; it increases more rapidly for groups than for rich 
clusters because of shallower gas-density slope in groups. 
We add the observed stellar fraction and the diffuse in- 
tracluster light to the gas fraction to obtain the total 
baryon fraction. We find that the average baryon fraction 
for all groups and clusters with M500 > 5 x lO^'^/ifj^AfQ 
increases steadily with radius, reaching the cosmic value 
and becoming flat as a function of mass when measured 
within the virial radius (for the LCDM cosmology). This 
suggests that baryons are not missing in clusters, but are 
simply located in cluster outskirts. This picture is con- 
sistent with heating processes (such as shock-heating of 
the intracluster gas, as well as supernovae and AGN feed- 
back) causing the gas to expand to the cluster outskirts. 
Upcoming observations in the X-rays and SZ should be 
able to detect these baryons. 
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FIG. 1: The observed cluster gas and baryon fraction within 
i?50o as a function of cluster mass, Msoo- The observed cos- 
mic baryon fraction is shown by the shaded band (l-cr; [14|V 
(a) Observed cluster eas fraction from Chandra and XMM- 
Newton within 7?5oo [a-TlH. The average for the four mass 
bins is shown by the filled circles, (b) Observed cluster gas 
and baryon fraction within -R500 for the averages of the four 
mass bins. Error bars are l-cr errors on the mean. Horizontal 
bars represent the mass range of each bin. 
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FIG. 2: Average observed gas density slope jSgas (where 
pgas oc r"^'''"'=) at > i?500 from ROSAT, Chandra, XMM- 
Newton, and Suzaku [3, l3ll - [34| as a function of cluster tem- 
perature. Filled circles are bin-averages of 39 ROSAT clus- 
ters [sJl; filled diamonds are ROSAT averages of hundreds 
of stacked optical clusters [32 . their richness classes and 
1]; empty circles are bin-averages of 10 Chandra clusters [g|; 
empty diamond is Suzaku observation for Abell 1795 [3^ : 
empty square is XMM observation for A1983 ^3f|]. The dashed 
line illustrates the average trend of Pgas with T. Error bars 
are 1-cr errors on the mean. Horizontal bars represent the 
temperature bin size. 
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FIG. 3: Cluster gas and baryon fraction (3a and 3b, respec- 
tively) within _R5oo (observed), _R200 and Rvir{~ Rioo) (ex- 
trapolated using observational data as discussed in the paper). 
The gas and baryon fractions are presented as a function of 
cluster mass for the four averaged binned samples. The error- 
bars are the l-cr error on the mean for each bin; horizontal 
bars represent the mass range of each bin. The observed cos- 
mic baryon fraction is presented by the shaded band (1-a). 
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FIG. 4: The dependence of the baryon fraction on radius, 
from -Rsoo to Rvir{= Rioo), for two representative mass bins 
(bins 2 and 4; the others show a similar trend). The radius 
is presented in units of -R200. The slow but steady increase 
of the baryon fraction with radius is apparent, reaching the 
cosmic value near R^ir for all cluster masses. The error-bars 
are 1-cr errors on the mean. 



